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System for Encoding and Manipulating Models of Objects 
Field of the Invention 

The present invention relates generally to computerized systems that 
generate models from the geometry, shape, and other attributes of objects, 
5 and, more particularly, to systems that manipulate the models. The models 
are usable by computer-aided design tools, robotics, molecular modeling, 
kinematics and dynamics of Unked bodies, colHsion detection, satellite 
positioning, articulated characters in animation and for use in the 
educational, scientific, business, and entertainment fields. 

10 

Background 

Points, lines, planes, and spheres of all dimensions can be used to describe 
basic objects or components. These basic objects can be assembles into more 
15 complex objects. Therefore, the manipulation and proximity testing of the 
basic objects has wide-ranging implications for many problems in 
engineering, computer graphics, and the practical sciences. 

For example, in robotics and many haptic devices, a tool is typically 
20 composed of rigid links connected to each other at joints. The possible 
positions of the distal end of each link, assuming complete rotational 
freedom at the joints, is of course a sphere, see U.S. Patents 6,070,109, 
6,064,168, and 5,973,678. 
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In order to check for interference with an obstacle in the vicinity of the tool, 
one can first check to see if bounding spheres intersects. This is a much 
simpler computation than determining exactly the coordinates of all of the 
individual parts. Only when the bounding spheres intersect does it become 
5 necessary to do the more rigorous computation on the exact coordinates to 
determine actual interference. 

In rigid body kinematics and CAD/CAM, a frequent problem is to check for 
collisions between arbitrarily shaped objects, see U.S. Patents 6,054,997 and 
10 5,943,065. This problem can also be simplified by first circumscribing each 
object with a bounding sphere. Now, one can first simply check for 
collisions between spheres, and only when there is a possible collision of 
spheres does one have to do the more complex computation to check for 
actual collision. 

15 

In graphics rendering, it is common to parametrically represent objects as 
polygons or Bezier surface patches. Frequently, it is necessary to rotate 
objects. The rotation of a complex object can simply be done by rotating 
points on the patches (planes) along spherical surfaces having their origin at 
20 the center of rotation. 

In fluid dynamics, and electrostatic problems, one frequently manipulates 
angle preserving conformal surfaces, see U.S. Patent 5,453,934. In higher 
dimensional geometries, a conformal surface can be represented by spherical 
25 rotation. Clearly, operations with a spherical surface would be much easier 
than operating on a complex surface such an airfoil. 
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In a Cartesian model of Euclidean space, the calculations performed on 
spheres are generally of the form: 

Operator((x^hlf +(y+ilf+(z-^jl)-rl^, (x+hlf +{y+i2f+{z^j2)-r2^) 
where "Operatof represents some spherical operation such as intersection, 
5 union, containment, tangent, translation, rotation, etc., and combinations 
thereof. In most of these application, event these simpler spherical 
computations, due their enormous numbers, still consume most of the 
system's resources when manipulating object models. Similar calculations 
can be performed on planes. 

10 

Therefore, it is desired to provide a new method for representing models of 
objects which solves these problems for practical modeling and simulation 
problems of physical objects. Furthermore, it desired to manipulate the 
models without having to consider the actual Cartesian coordinates of the 
15 underlying objects. 

Summary of the Invention 

The present invention encodes points of components of objects measured in 
20 Euclidean three-space G(3,0) to null vectors in a general homogeneous four- 
space G(4,l). This encoding makes it possible to directly determine scalar 
distances between points by vector inner products. 

In addition, planes and spheres planes 'mG(3) can be represented by non- 
25 null vectors in G(4J} so that distances between planes can also be computed 
by vector inner products which yields a scalar. 
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All relations among points, lines, planes, and spheres inG(3) can be 
represented or determined in G(4J) without having to resort to their actual 
coordinates, as in the prior art. 

5 Common operations such as 3-D rotational and translational kinematics of 
objects is described by a single spinor equation in G(4J ). The 3-D rotational 
and translational dynamics of a rigid object is described by a single bi vector 
equation in G(4J). 

10 More particularly, the invention provides a method for modeling an object 
composed of one or more components. Data are input to a memory of a 
computer system for each component of the object. The data of each 
component point include Cartesian coordinates expressed in Euclidean space 
of a plurality of points x of each component. 

15 

Each component point x is encoded as a vector x in a general homogeneous 
space hyx = (x+2X^e + e^)E = xE-2^e + , where e and e^^ are basis 
null vectors of a Minkowski space E. General homogeneous operators are 
associated with each data point to generate a model of the object. 

20 

The general homogeneous operators are applied to each encoded point of the 
associated component for each component to manipulate the model of the 
object. The operators can be in terms of processes executing in the computer 
system. 
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Brief Description of the Drawings 

Figure 1 is a block diagram of a modeling system according to the invention; 

5 

Figure 2 is a graph of a Minkowski plane used by the present invention; 

Figure 3 is a graph of a homogeneous space where objects are modeled 
according to the invention; 

10 

Figure 4 is a table of transformation operations to be applied to objects 
according to the invention; and 

Figure 5 is a schematic of a robot and an obstacle modeled by the invention. 

15 

Detailed Description of the Preferred Embodiment 
System Overview 

20 As shown in Figure 1, our invention is a system 100 for generating and 
manipulating a computerized model of a real or virtual object 101, for 
example a robot. The object is composed of basic components such as links 
and joints. The object is the "input" 102 into the system 100. The input can 
be described in terms of the geometry (data or constructs) 103 of the basic 

25 components, and the operations (methods or operators) 104 that can be 
applied to the basic components. 
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The data 103 can be points, lines, planes, and spheres. The data can be 
acquired by measuring the objects, or the data can be parametrically defined. 
The data are stored in a memory of a computer system. The methods 104 can 
be linkage dynamics, torque forces, transformations, rotations, constraints, 
5 and the like. The methods can be implemented as programs executing in the 
computer. 

The system 100 according to our invention operates the model in a 
homogeneous space (G(4J)). Homogeneous space is described in greater 

10 detail below. We presume that the methods 104 are already in a 

homogeneous form. However, the data 103 expressing the geometry of the 
basic components most likely use Cartesian coordinates of Euclidean space 
{G(3,0)). Therefore, step 115 can convert the data 103 to homogeneous form 
105 where necessary. The encoder 110 combines the homogeneous data and 

15 corresponding homogeneous methods to generate classes that together for an 
object-oriented programming structure (OOPS) 130. 

The modeler 120 operates on the OOPS 130 using run- time parameters 121 
to determine distances, intersections, and tangencies of the basic 
20 components, and to perform operations such as rotations and displacements 
of the basic components. In an object oriented programming structure, 
methods are applied to data. Output of the modeler 120 can be rendered as 
images on a display device 109. 

25 As will be described, the result of many of the operations, for example, 

operation that measure distances, result in a scalar. Therefore, the result does 
not need to be decoded, and can be directly used in the context of the 
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originating Euclidean space. In other words, as an advantage of our 
invention, the output 108 of the modeler 120 can be directly processed by 
any standard rendering engine, or the robotic tool itself, without decoding or 
converting. 

The homogeneous constructs and operators (OOPS) 130 generated by the 
encoder 110 and used by modeler are derived using geometric algebra. 
Geometric algebra is described extensively by Hestenes in "New 
Foundations for Classical Mechanics'' Kluwer (second edition), 1999. 
Geometric algebra has been used extensively to describe the mechanics of 
particles, rigid bodies, biomechanics, robotics, computer vision, computer 
aided design, orbital mechanics of space programs, and many other practical 
fields. Geometric algebra is a general modeling system for a wide range of 
physical phenomena. Starting from a small number of simple rules that 
define the geometric algebra, one can derive vector algebra, quatemions, and 
spinor representations of quantum mechanics. However, applications of 
geometric algebra in our "generalized" homogeneous space are new. 

The standard model for Euclidean space E"^ is an n-dimensional real vector 
(Cartesian) space or, equivalently, a set of real coordinates. One trouble 
with this model is that the origin is a distinguished element, whereas all the 
other points of are identical. This deficiency in the vector space model 
can be corrected by removing the origin from the plane and placing it in a 
space one dimension higher. 
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We do this by placing our model in general homogeneous space. In general 
homogeneous space there are no coordinates. Our model also provides for 
an adequate representation for Euclidean points or lines at infinity. 

The idea underlying our general homogeneous space for "points" in is to 
remove the troublesome origin by embedding i?'' in a space of higher 
dimension with a displacing "null" vector. As stated above, we measure or 
model objects in the Euclidean space E"" as a set of points with parametric 
properties. A standard way to do this is by identifying each Euclidean point 
of an object with a vector x in R^, as expressed by the isomorphism 

We call this the inhomogeneous model of because its origin Oisa 
distinguished point in even though all points in are supposed to be 
identical. 

Constructs 

Following are some of the basic data constructs used by our encoder and 
modeler. 

Scalars and Vectors 

Geometric algebra is built over a real, n-dimensional vector space. Thus, we 
use the structures of scalars and vectors as are well known in the art of 
geometric algebra. 
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Null Vector 

If a vector a^O, but the absolute value of vector \a\ = 0, then vector a is a 
null vector. In geometric algebra a null vector can exist when a^O. 

5 

K-blade 

The outer product of k vectors, determined by 
a]Aa2A . . . Aajc^ 

10 s defined as the antisymmetric part of the geometric product aja2 . . , Uk. The 
resulting quantity is called a k-blade, where k is called the grade of the 
blade. 

Pseudoscalar 

15 

The largest grade k-blades are called pseudoscalars. We reserve / as a 
special symbol for a largest k-blades because these are all scalar multiples of 
a single unit pseudoscalar. 

20 Versor 

A quantity 

U=UkUk_j . . ,u (1) 
generated by a product of non-null vectors is called a versor 

25 
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Operators 

In order to properly define a class, it is also necessary to define the 
applicable operators (methods) 104. For the above constructs, the operations 
are defined according to geometric algebra. Thus, for example, a linear 
transformation U on the vector space i?^'^ is said to be orthogonal when the 
transformation preserves the inner product of vectors, as specified by: 
(UaHUb) = a^b. 

Every orthogonal transformation U is determined by a versor U via the 
equation: 

U{x) = zUxir^, (2) 
where 8 is the parity of C/. This reduces orthogonal transformations to 
multiplication by versors. 

A multiplication by versors has at least two great advantages. First, it 
reduces composition of orthogonal transformations to multiplication of 
versors. Second, it enables us to represent and visualize any orthogonal 
transformation as a set of vectors via equation (1). Third, the result is a 
scalar that provides the answer to the problem originally expressed in 
Euclidean space. 

In addition, we also define a new class, which leads to the encoding of the 
input 102 into new data structures 104 that are subsequently operated on 
with new methods. 
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Minkowski plane E 

As shown in Figure 2, a Minkowski E plane 20 lis generated, from two 
unique basis null vectors e 202 and e^^ 203, i.e., a basis vector {e,e,} with a 

5 fixed direction and variable scale. The null basis vectors have the following 
constraints 

e =0 = e.i. , 

and 

= 1 + e ^£ ~ I - E, 

10 

where E = eAe^, is a unit 2-blade, that is, = 1. This is called the 
Minkowski plane E. An alternative but equivalent orthonormal unique basis 
vector, with a variable direction and fixed scale, is {e.,e+} 204. Because the 
orthonormal basis vector 204 has a fixed scale but an arbitrary direction, we 
15 prefer the null basis vectors 202-203. 

Conformal Split 

As shown in Figure 3, a conformal split allows us to relate objects and 
20 operations between spaces of different dimensions, particularly to relate 
operations (methods) on data expressed in Euclidean and homogeneous 
spaces. The conformal split is realized by intersecting the cone of null 
vectors 301 generated by {e,e,} with the plane orthogonal to e, resulting in 
the so-called horosphere 300. The horosphere was originally described by F. 
25 A. Wachter (1792-1817), a student of Gauss. He showed that a certain type 
of surface in hyperbolic geometry is geometrically equivalent to Euclidean 
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Space. Our invention shows that the horosphere 300 constitutes a 
homogeneous model of Euchdean space. 

Unhke Wachter, we formulate the horosphere in terms of geometric algebra. 
Thus, our horosphere finally becomes a practical tool many real-world 
modeling applications. 

Homogeneous Point 

Using the horosphere 300, we encode any input data point x of an object 
component in the data 103 in standard Euclidean space, into a new point x in 
a homogeneous space, using 

X = (x + 2X^^ 4- e^^E = xE- 2X^e + e.-, , (3) 
which expresses x in terms of input x. 

Homogeneous Point, Line, Plane and Hyperplane 

For independent input vectors ao, a],. . . a^ of the data 103, where m=0 
defines a point, m=l defines a line, m-2 defines a plane, m=3 defines a 
hyperplane, etc., a corresponding vector p is encoded respectively for a 
homogeneous point, line, plane, hyperplane by 
p = IE (cA aoA a] A... Aam). 
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Homogeneous Sphere 

An input sphere of a radius r, and a center c iiaving coordinates of m 
dimensions, we encode a m-dimensional vector s, called the homogeneous 
sphere, by 

C + 27 e. 

Note that in our new homogeneous classes all objects, i.e., points, lines, 
planes, hyperplanes, and spheres of any dimensions, are exclusively 
represented as vectors. 

Thus, our invention leads to much simpler novel operations used by our 
modeler 120. Specifically we describe the following new methods, which 
operate on our homogeneous representation of the basic homogeneous 
constructs. 

Distance between Points 

The Euclidean distance dab between two Euclidean input points a and b is 
determined for homogeneous encoded points a and b by 

2 2 

dab =(a-b) = -2a»b, 
a scalar directly useable in Euclidean space without having to do any further 
conversion. 
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Length of a Line Segment 

The length of line, or line segment through two input points a and b is 
determined for encoded homogeneous points a and b by 
(length)^ (eAaAbf = (a- bf^ 

also a scaler. 
Triangular Area 

The homogeneous area of plane given by input points a, b, c, or the 
triangular plane segment with vertices a, b, c is determined by 

(area)^ = 4 {eAaAbAcf^ 

a scalar. 

Distances between Points and Objects 

The distances between homogeneous point a, plane p and sphere s are 
determined by the respective inner products 
a^py a^b, a*s, p*s, 

a scalar. 

Distances between Spheres 

The scalar distances between two homogeneous spheres S] ^ Cj + 2T1 e and 
S2~ C2+ 2r2 eis determined by 

Si^S2 = C1-C2 + 2 (r;^ -h r2^) = -2 {{ci ~ C2f - (rj^ + r2^)]. 
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Rigid Body Motion 

According to Charles' Theorem, any rigid displacement can be expressed as 
a screw displacement. Therefore, the motion of a rigid body in homogeneous 
space is determined by a time dependent displacement versor D = D{t) 
satisfying the differential equation 

b =iyA 

with a "screw velocity V given by 

where co is the angular velocity and v is the translational velocity of the 
body. 

Dynamics of a Rigid Body 

The dynamics of a rigid body in homogeneous space is determined by the 
differential equation 
P = W 

with 

P = -/L + e.p and W=-/T + e.F, 
where L is the angular momentum and p is the translational momentum of 
the body, while T is the net torque and F is the net force on the body. 

Conformal Transformations 

When equation (3) is substituted into (2) we get 

8[/(x +\x^e + e:)Eir'E = o(x' + + ^.,) 
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where a is a scale factor. The scale factor can be used for normalization 
purposes. This equation determines a transformation of each input vector x 
into a unique vector 
x' = g(x). 

This induced function g(x) is a conformal transformation of input vectors 
and every conformal transformation can be generated in this way. A 
transformation is said to be conformal if it preserves angles between 
intersecting curves in the space being transformed. 

Every conformal transformation function g(x) can be represented by a versor 
U and composition of conformal transformations is reduced to multiplication 
of versors as expressed by equation (2). 

The ability to compute conformal transformations by multiplication this way 
is new and important innovation for manipulation of objects. 

Figure 4 shows all conformal transformations that can be applied to input 
points and homogeneous points. In our homogeneous form, composition of 
transformations is accomplished by simple linear multiplications, as opposed 
to nesting functions in the standard Euclidean form. Also, note that our 
homogeneous transformations never require division or addition. 

In Figure 4, column 401 shows different types of transformations, column 
402 non-linear transformations of prior art EucHdean space, column 403 
linear homogeneous transformations, and column 404 scalars. In Figure 4, 
the symbols are used as follows, x a point in Euclidean space, n normal to a 
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plane, e a null vector, ^5 a displacement, p a radius, c a center, R a rotation, T 
2i translation, a and X. 

Application Example 

Figure 5 shows a robot 501 operating within the confines of an obstacle 502. 
The obstacle is associated with bounding planes and box 503, and bounding 
spheres 504, and the robot is associated with bounding spheres 505. As the 
robot 501 moves, it is important to avoid hitting the obstacle 502. The 
calculation to determine interference must be made quickly. 

Our invention models the robot and obstacle by a few well chosen spheres 
and planes in homogeneous form. We then test for interference using the 
homogeneous operations as defined above. If the spheres that model the 
robot are not allowed to intersect the spheres and planes of the obstacle, then 
it is ensured that the robot will not collide with the obstacle. Therefore, we 
encode the robot and its operation in homogeneous form, and then apply the 
appropriate operations (methods) 104 to determine the distances between the 
various spheres, planes, and lines. 

In addition, our operations can determine the dynamics of each robot link, 
for given input torques, forces and momentums, and link lengths. Similar 
methods apply, for instance, in CAD/CAM where fast intersection of objects 
can use a collection of bounding spheres. 

Finally, the motion (transformations) of the robot can be determined because 
all rigid body motions are conformal. The modeler 120 can determine the 
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motion by a linear sequence of multiplications, instead of the more 
traditional combinations of adds and multiplies. 

Conclusions and Effect of the Invention 

We have described how an object, measured or defined in Euclidean space, 
can be encoded as a model in a homogeneous space. After encoding, 
complex non-linear operations in Euclidean space become simple linear 
operations in homogeneous space.. This greatly simplifies operations on 
planes and spheres that are frequently used in robotics, kinematics, and 
computer graphics. 

Although the invention has been described by way of examples of preferred 
embodiments, it is to be understood that various other adaptations and 
modifications may be made within the spirit and scope of the invention. 
Therefore, it is the object of the appended claims to cover all such variations 
and modifications as come within the true spirit and scope of the invention. 
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Claims 

1. A method for modeling an object composed of one or more components, 
comprising: 

inputting data for each component of the object, the data including 
Cartesian coordinates expressed in Euclidean space of a plurality of points x 
of each component; 

encoding each point x as a vector x in a general homogeneous space 

1 2 1 2 

byx = (x+2X^ + e.)E = xE - 2^ e + e,-, , where e and e^^ are basis null 
vectors of a Minkowski space E; and 

associating a plurality of general homogeneous operators with each 
data construct to generate a model of the object. 

2. The method of claim 1 further comprising: 

supplying run-time parameters for the plurality of operators; and 
applying the plurality of general homogeneous operators to each 

encoded point x of each associated component to manipulate the model of 

the object, 

3. The method of claim 1 further comprising: 

measuring a scalar dab between two component points a and b 
encoded as general homogeneous points a and b by dab^ = (a~ bf = -la^b. 
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4. The method of claim 1 wherein a Hne through component points a and b 
encoded as general homogeneous points a and b is modeled by eAUAb, and a 
length lab of a line segment connecting component points a and b is 

(i^^)^ = {BAaAbf -{a- bf. 

5. The method of claim 1 wherein a plane through component points a, b, 
and c encoded as general homogeneous points a, b, and c is modeled by 
eAUAbAc, and an area Aabc defined by component points a, b, and c is 

(Aabc) = 4 (eAaAbAc) . 

6. The method of claim 1 wherein a sphere s with radius r centered at a 
component point c encoded as a general homogenous radius r and center c 

is generated by a vector s = c + 2^ e. 

7. The method of claim 1 wherein a sphere s determined by four component 
points a, b^ d encoded as general homogeneous points a,b,c,d is generated 
by s = /E {uAbACAd), where / is a largest k-blade. 

8. The method of claim 7 wherein one of the general homogeneous points 
a,b,c,d is equal to the point e so that s defines a plane through the point e. 

9. The method of claim 5 wherein a distances between a component point a 
and a component plane p is an inner products a^p of an encoded point a and 
an encoded plane p. 
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1 10. The method of claim 6 wherein a distances between a component point a 

2 and a component sphere s is an inner product a*s of and encoded point a and 

3 the encoded sphere p. 

1 11. The method of claim 6 wherein a distance between two component 

2 spheres Si and S2 encoded as spheres Sj = Cj 2 rj^e and ^2 = + 2 is 

3 generated by srs2 = crc2 + 2 (^/^ + r2^) ^ ~A [(c; » ^2)^ ™ (r/ + r2^)]. 

1 12. The method of claim 1 wherein the object is a rigid body, and a motion 

2 of the rigid body is determined by a time dependent displacement versor 

3 D-D(t) satisfying a differential equation/) = 2 VD, with "screw velocity" V 

4 given by y = -/CO + ev, where 00 is a velocity and v is a translational velocity 

5 of the rigid body. 

1 13. The method of claim 12 wherein dynamics of the rigid body are 

2 determined by a differential equation P -W, where P - ~IL + e^.^, and 

3 W= -/T + ej^, where L is an angular momentum and p is a translational 

4 momentum of the rigid body, while T is the a torque and F is a net force on 

5 the rigid body. 

1 14. The methods of claim 12 wherein the rigid body includes n linked rigid 

2 components, and a motion of the rigid body is modeled by n time dependent 

3 displacement versors Du Z)2, • • - , , and a motion of a k^^ hnked rigid 

4 component is determined by a versor product D1D2 . . . Dk- 
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1 15. The method of claim 1 wherein the objects is a robot composed of a 

2 plurality of links connected at joints. 
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Abstract of the Disclosure 

A method models an object composed of one or more components. Data are 
input to a memory of a computer system for each component of the object. 
The data include Cartesian coordinates expressed in Euclidean space of a 
plurality of points x of each component. Each component point x is encoded 

as a vector x in a general homogeneous space by x = (x + + e.^E -xE~ 

1 2 

2X e + e^i., where e and e-^ are basis null vectors of a Minkowski space E. 
General homogeneous operators are associated with each data point to 
generate a model of the object. The general homogeneous operators are 
applied to each encoded point of the associated component for each 
component to manipulate the model of the object. 
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DECLARATION AND POWER OF ATTORNEY 

As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated 
below next to my name , 

I believe, the below named inventor is the original, first 
inventor of the subject matter which is claimed and for which a 
patent is sought on the invention for System for Encoding and 
Manipulate Models of Objects, the specification of which is 
attached. 

I hereby state that I have reviewed and understand the contents 
of the above-identified specification, including the claims. I 
acknowledge the duty to disclose information which is material to 
patentability in accordance with Title 37, Code of Federal 
Regulations , §1.56. 

I hereby declare that all statements made of my own knowledge are 
true and that all statements made on information and belief are 
believed to be true; and further that these statements were made 
with the knowledge that willful false statements and the like so 
made are punishable by fine or imprisonment, or both, under 
Section 1001 of Title 18 of the United States Code and that such 
willful false statements may jeopardize the validity of the 
application or any patent issued thereon. 

1 hereby appoint the following attorney with full power of 
substitution to act exclusively for me to prosecute this 
application and transact all business in the Patent and Trademark 
Office connected therewith: Dirk Brinkman, Reg. No. 35,460. 

All correspondence should be addressed to: 
Dirk Brinkman Esq. 
Patent Attorney 
2 01 Broadway 

Cambridge, Massachusetts 0213 9 

All telephone calls should be directed to Dirk Brinkman, 
telephone number (617) 621-7539. 
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Inventor's Full Name: 


Alyn Rockwood 


Inventor's Signature: 




Date : 




Residence : 

(City, State and/or 
country) 


18375 Highland Estate DR 
Colorado Spring 
Colorado, 80908 


Citizenship : 


United States of America 


Post Office Address: 


18375 Highland Estate DR 
Colorado Spring 
Colorado, 80908 



Inventor's Full Name: 


Hongbo Li . My current address : 


Inventor ' s Signature : 




Date : 




Residence : 

(City, State and/or 
country) 


Building 902, Room 1002 
Zhong Guan Cun District 
Beijing 100080 
P. R. China 


Citizenship : 


P.R China 


Post Office Address: 


Building 902, Room 1002 
Zhong Guan Cun District 
Beijing 100080 
P, R. China 



Inventor's Full Name: 


David Hestenes 


Inventor's Signature: 




Date: 




Residence : 

(City, State and/or 
country) 


2416 S. Palm Drive 
Tempe, Arizona 85282 


Citizenship : 


United States of America 


Post Office Address: 


2 416 S. Palm Drive 
Tempe, Arizona 852 82 



2 



